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Abstract 



>0 I Criticality in a fluid of dielectric constant D that exhibits Ising-type be- 

pg I havior is studied as additional electrostatic (i.e., ionic) interactions are turned 

o, 

0^ I on. An exploratory perturbative calculation is performed for small ionicity 



as measured by the ratio of the electrostatic energy e^/D a (of two univa- 
lent charges, ite, separated by the atomic/ionic diameter a) to ks T^ which 
C3 ' represents the strength of the short-range nonionic (i.e., van der Waals) in- 

o' 

^ ' teractions in the uncharged fluid. With the aid of distinct transformations 

K^ \ for the short-range and for the Coulombic interactions, an effective Hamilto- 

H 

C^ I nian with coefficients depending on the ionicity is derived at the Debye-Hiickel 

limiting-law level for a fully symmetric model. The crossover between classical 

(mean-field) and Ising behavior is then estimated using a Ginzburg criterion. 
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This indicates that the reduced crossover temperature depends only weakly on 

the ionicity (and on the range of the nonionic potentials); however, the trends 

do correlate with the, much stronger, dependence observed experimentally. 
61.20.Qg, 05.20.-y, 05.70.Jk 
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I. INTRODUCTION 

Early studies of the liquid- vapor critical point of ionic systems were carried out on molten 
salts and related systemsJlTQ with relatively poor temperature stability and resolutionD 
Results for ionic critical behavior at low reduced temperatures were first published by Pitzer 
and co-workers a decade ago.l3 Subsequently other solutions of organic salts, with low critical 
temperatures, were studied yielding somewhat puzzling results. 

One of the systems studied in detail was triethyl-n-hexylammonium triethyl-ra- 
hexylboride (N2226^B2226) in diphenyl ether. In 1990, Singh and PitzerQ reported a 
coexistence-curve exponent, /? ~ 0.476, close to the mean field value of |. Their data 
are compatible with Ising criticality (/? ~ |) only if the mean field to Ising crossover occurs 
at a reduced temperature, tx = |T — Td/Tc ~ 10~^.I§Q This value of tx is unusually small — in 
ordinary fluids tx ~ 1 and the observed critical exponents are clearly of Ising character. In 
1992, Zhang et alB confirmed the mean- field critical behavior of this system by reporting a 
susceptibility exponent, 7 ~ 1.0. Q 

A different system also studied nearly a decade ago was tetra-n-pentylammonium bro- 
mide dissolved in water. Japas and Levelt SengerstO reported an Ising-type coexistence curve 
with j3 = 0.319-0.337. Other studies reported Ising or mean-field criticality depending on 
the particular ionic solution. For reviews of the various experiments, see: Refs.u'0Mll3 

The presence of the solvent, which lowers the critical temperature of the salts, complicates 
the situation considerably. The bare Coulomb interaction is reduced by a factor that depends 
on, D, the static dielectric constant of the solvent. Additionally, in polar solvents, solvation 
shells are formed around each ion and the resulting effective interactions between the ions are 
not well characterized at present. The mechanisms driving criticality in these various systems 
are not clear, but there appear to be two classes of critical behaviour: Coulombic systems, 
where the electrostatic interactions primarily drive the phase transition and solvophobic 
mixtures in which the phase transition is mainly driven by the short-ranged interactions 
that favor solvent-solvent and solute-solute contacts.Q'Qii^Jl^ 



In Table | we display data on the melting temperature of three salts and their critical 
temperatures in solution. Clearly, if T^ is lower than Tc the Coulomb interactions are 
too weak to drive the phase transition, especially in high dielectric-constant solvents.t3 
However, there are cases where the ionic interactions are strong enough to play a role at 
temperatures of the order of T^. In these systems,t3 a crossover from mean-field to Ising 
criticality was observed, with reduced crossover temperatures tx , that decrease as the relative 
ionic character or ionicity increases, i.e., the dielectric constant of the solvent decreases. Ising 
exponents are observed only closer to Tc as the ionicity increases. It is plausible to suppose 
that the effect of the Coulomb interactions is the reduction of the crossover temperature; 
this reduction might be so strong that Ising critical behavior could not be observed in 
practice. Ideal Coulombic electrolyte criticality would then be characterized by classical 
(mean-field) behavior, while solvophobic electrolyte criticality remained in the 3D Ising 
universality class.cJa 

However, a theoretical description of the dependence of tx on the Coulombic and other 
interactions is still lacking. Indeed, previous attempts to study the crossover temperature of 
model ionic fluids do not seem to predict any reduction, at least for the restrictive primitive 
model (RPM) consisting of an equal number of equisized positive and negative hard-spheres. 
Following the sugestion of Fisher and Levin,t3 Leote de Carvalho and Evansli3 and Fisher and 
Leet3 employed a Ginzburg criterion which led them to conclude that the reduced crossover 
temperature of the RPM is similar to that of ordinary fluids. 

Although the study of criticality in the RPM has proceeded for half a decadeQcjEZl 
progress has been slow and no definitive conclusions have yet been reached. Fisher and 
LevinllB extended the original Debye-Hiickel theory and established the importance of in- 
cluding dipolar Bjerrum pairs^^ and of accounting for the solvation of the neutral, but elec- 
trostatically active, pairs in the screening ionic fluid. These authors have obtained the best 
theoretical estimates of the critical density and temperature of the RPM when compared 
with the results of computer simulations.c3E3iH3 

However, the direct Fisher-Levin approach does not address the question of the univer- 



sality class of the RPM. On the other hand, computer simulations of the critical behavior 
of the RPM have been carried out; but these also have not provided clear-cut answers. Cal- 
liol et alB^ have obtained evidence of consistency with Ising exponents, as have Orkoulas 
and Panagiotopoulos;c3'E3 but Valleau and TorriecH claim that no sign of Ising behavior is 
observed in their simulations above criticality. 

Nevertheless, it is instructive to pose a more modest question, namely: "How does the 
introduction of weak electrostatic interactions affect the near-critical behavior of a fluid 
that, when uncharged, displays Ising criticality ?" While one cannot discard the possibility 
that the long-range Coulomb interactions discontinuously change the universality class,Q the 
Debye screening of the ions at any non-zero ionic strength can be regarded as giving rise 
to an effective short-ranged interaction,LH13Ey and thus — at least in a perturbative sense 
— we may suppose that weak electrostatic couplings will not change the character of the 
transition. They must, nonetheless, affect all the non-universal cnticaA properties: and many 
of these are of independent interest. 

Accordingly, in this article, we address this last issue. Specifically, we consider an ordi- 
nary fluid of dielectric constant D with only short-range, or van der Waals interactions, that 
exhibits criticality of Ising character at a temperature T°. For simplicity in this essentially 
exploratory study, we focus on a fully symmetric model in which the particles of the original, 
uncharged fluid are all identical having repulsive cores of diameter a and fixed attractive 
pair interactions, the strength of which is measured by /c^T^. 

Then we ask how the reduced crossover temperature, tx, changes when electrostatic 
interactions are switched on so that half the particles become positive ions of charge +e 
and the rest negative ions carrying charge — e. The strength of the Coulomb forces may 
be measured by the "contact energy" e'^/D a, so that one can imagine either increasing e 
from zero or, for fixed charges, decreasing the dielectric constant D from very high values. 
Naturally, the critical temperature of the system will change and it is thus reasonable to 



define the ionicity of tlie critical cliarged system by 

e^ 1 

X= - = _ (I) 

wfiere T* = D aks T/e^ is tlie standard definition of reduced temperature in primitive model 
electrolytes.Q'0 The ionicity will constitute our perturbative parameter. 

To treat the model, we first consider, in Section II, the purely ionic limit (where kB T^ = 
0) and, using a lattice cutoff, reformulate the Coulomb interactions via a quadratically trun- 
cated sine-Gordon transformation.c3'E3 This corresponds to a Debye-Hiickel limiting-law level 
of analysis although the expansion could be carried to higher order.B^ Then, in Section III, 
we introduce (charge independent) short-range interactions with the aid of a Kac-Hubbard- 
Stratonovich (KHS)e3'EI transformation. In the resulting (approximate) Hamiltonian we 
integrate out the sine-Gordon charge field iP{t), perturbatively to obtain an effective Landau- 
Ginzburg- Wilson (LGW) Hamiltonian in terms of a spin (or density) field, s(r), but with 
coefficients depending on the ionicity, I. Finally, effectively truncating at the s^ level, we 
estimate the crossover temperature with the aid of a Ginzburg criterionli3li3E3 in Section IV. 
This treatment indicates that, in general, the reduced crossover temperature decreases by 
an amount of order 25%. Such a weak change is far too small to describe the experimental 
observations although the overall trends seem to correlate: see Table II and the discussion 
and conclusions in Section V. 

II. LATTICE CUTOFF MODEL FOR THE IONIC FLUID 

Consider a system of unit-valence positive and negative ions with repulsive cores but 
otherwise interacting only through Coulomb forces. To account approximately for the ionic 
cores we may suppose the ions reside on the sites of a lattice of spacing a equal to the atomic 
diameter. Alternatively, in momentum space, we can employ a cutoff on a corresponding 
Brillouin zone or, for convenience, approximate this by a sphere of radius vr/a. The grand 



partition function for the lattice system with sites labelled i, j, etc., is 

Z = y^ exp< — /i+ n_|_ — /i_ n_ ^ x 

r 1 1 ^^^ 

{p+,p-} «>i 

where the first sum runs over ra+ and n_, the total numbers of positive and negative ions, 
subject to 72+ + r7,_ < A^, the number of lattice sites, while the second sum runs over different 
lattice configurations, for given numbers, r2+ and n_. The occupation variable p+{i), is 1 
if a positive ion is at site i, but otherwise, with a similar definition of p_. We omit the 
possibility of ionic pairing (molecule formation)!] so that a lattice site cannot be occupied by 
more than a single ion. The remaining symbols have their usual meanings: /i+ and fi- being 
the chemical potentials (divided by ks T) for the positive and negative ions, respectively. 

Since this is a first approach and we are interested in real fiuids (rather than details 
induced by the lattice structure) we do not take the Coulomb interaction proportional to 
the lattice Green's function but simply use the reduced three-dimensional isotropic Coulomb 
potential so that in @) we have 

where Vij is the distance separating sites i and j. 
Given a positive definite matrix Aij, the identity 



T1 V^'l l-l 1 






(where i = v^— T) may be used to rewrite the grand partition function of a repulsive system 
with discrete variables in terms of the partition function of a system with continuous vari- 
ables: this constitutes the well known sine-Gordon transformationEa that may be regarded 
as a modified form of the Kac-Hubbard-Stratonovich transformation used for attractive 
interactions.Ey'El These transformations are standard and need not be discussed in further 
detail here. 



The diagonal elements of the electrostatic interaction matrix, (fdhi), are chosen to 
guarantee that (pc{i,j) is positive definite,c^ but are otherwise arbitrary since a lattice site 
cannot be occupied by more than a single ion. Thus (^) may be used to write (|^) as, 

z = I ^ exp{-^-j2mvc\hj)^{j)}m, (5) 

where the weight function is given by 

Wo=X] Y^ exp|-i^[p+(j) -p_(j)]^(j) -/i+n+-^_n_|. (6) 

n+,n. {p+,p-} J 

The continuous field ipi^j), ranging from — oo to +00, is conjugate to the (reduced) discrete 
charge density 

qij)=pAj)-P-ij)- (7) 

On the lattice one has Dtp = Y[~ d'ipiJ)^ while Zq is an unimportant constant given by 



Defining the thermodynamic fields fi = (/i+ + /i-)/2 and 7 = (yU+ — yU_)/2, and noting 
that the total ionic density may be written as 

p{j) = P+ij) + P-ij)=Q'ij). (9) 

the weight function takes the more convenient form 

Wo= Yl exp|-i5^g(j)[^(j) + i7]-5^g'(j)/i}- (10) 

g=-i,o,i j j 

Carrying out the sum over q here, leads to 

Wo = n(l + exp{-i[V'(j) + 17] -P} + exp{+i[^(j) + 17] - /i}) 

(11) 
= exp< y^ ln[l + 2 exp(— /i) cos{ip{j) + 17) 



It may now be checked, using (j^) and (|Tl|), that global charge neutrality, namely (g) = 0, 
is obtained if 7 = 0, i.e., when the chemical potentials of the positive and negative ions are 
the same. From hereon we assume this to be so and set 7 = 0. 

Provided z = e~^ < i, we may expand the logarithm in the weight function Wo with 
respect to V'(j)- By truncating the expansion at second ordeJ^ we find 

(Z/Z^^)'/^^l + 2;. + ---, (12) 

where we have defined the Debye-Hiickel limiting-law partition function as, 

^"" ^ / ^ exp{-i Y.'^{^) {v~\i.2) + (r||^'^(^'^')] V^(J)} 

'^ 1 "' ^ (^^) 

= / ^ ^^P{~2 Xl^(^)^ii'^^'-^)^(-^)}' 

where S{i,j) is the Kronecker delta and '^LL^hJ) is the Debye-Hiickel limiting-law (DHLL) 
reduced effective potential as we will demonstrate. 

Now, Fourier transformation of the electrostatic potential @ on the lattice leads to 

Mk) = ^[l + 0{k'a% (14) 

where £b = e^/DksT = a/T* is the Bjerrum length (the distance at which two elemen- 
tary charges interact with an energy ksT). The 0{k'^a?) terms here arise from the lattice 
structure or, equivalently within our approximation, from the repulsive core of the ion-ion 
interactions. Since our subsequent analysis will utilize only the leading k ^ Q behavior, the 
effects of the lattice structure entering this way, via the Coulombic potential, should not be 
significant; thus we will neglect the 0{k'^a?) term in (0). Alternatively, one may regard the 
pure 1/fc^ behavior to be retained in (0) as defining (via inverse Fourier transformation) the 
Coulomb potential as viewed on the lattice in real space: this would mean that (fdhj) as 
stated in (|^) should be modified slightly on distance scales of order a. However, in a fuller, 
more precise treatment, of the lattice primitive model it is possible that such contributions 
could play a more crucial role;^ but that issue lies outside the scope of our present analysis. 



It follows from ([T3|) and (|T^) that the DHLL effective potential in Fourier space is given 
by 

47r£B(T) k'^ + An^BV 
where r] = 2z/{l + 2 z) with z = e~^ and, for convenience here and in most places below, we 
express k, is, etc., using length units for which a = 1. By differentiating the grand partition 
function (^ with respect to ^ and expanding the result in powers of ip, we find the total 
density of ions is given by 

Retaining only the first term in this series, which implies f] ^ p, and transforming (plSf ) back 
into real space, we obtain 

exp(-r/^z3) , . 

^LL{r) oc , (17) 



where C,d = l/V^Trf^P, is the standard Debye screening length,E3 thus demonstrating that 
the correct DHLL effective potential for ionic systems may be obtained from a quadratic 
expansion of the sine-Gordon weight function, in ip — a well known result .EaiH^ 

In the next section we derive, using similar techniques, a field-theoretic description of a 
system with additional short-range forces. As explained, we will focus on the situation where 
these forces drive the critical behavior. Then by integrating out the electrostatic interactions 
at the DHLL level we obtain corrections to the nonuniversal critical parameters, including 
an estimate of the crossover temperature, tx- 

III. SYMMETRIC IONIC MODEL WITH SHORT-RANGE ATTRACTIONS 

By adding to the lattice-cutoff model an attractive symmetric short-range pairwise po- 
tential identical for all ionic pairs, so that v'++(z,j) = (pt^{i,j) = '■p\^{i,i) = (ps{hj), the 
grand partition function becomes 

10 



(18) 



{q} i i,j 

where we have set 7 = (/i+ — ^_)/2 = 0, to ensure global charge neutrality. The sum on 
{q} runs over all possible configurations of ions on the lattice. 

The sine-Gordon transformation may be used as before on the electrostatic interactions 
in ([T8| ) to obtain a description in terms of the field, ip. In general, the short-range interaction 
matrix has positive and negative eigenvalues^ and thus the KHS transformation cannot be 
applied directly.Ey Here, we assume that —ips{i,j) is positive definite and use the identity 



to transform the short-range attractive interactions. After these transformations the grand 
partition function ([ISD reads 



Z= f ^^ exp{-i J2 h^) ^c\hj) Hj) + si^ ^:\^,J) sU)] } Wi, (20) 

where the field s{j) is conjugate to the total ionic density g^(j): see (^. The normalizing 
factor is now 

Z, . , P-'" (21) 

^det[y?c(«,j)]det[y?,(«,j)] 

and the new weight function is 

>Vi = exp|j]ln[l + 2exp[-s(j)-/^] cosV^(j)]}. (22) 

j 

Now let us shift (and rescale) the s field by putting 

Sj = l{sj + fi-ln2), (23) 



and expand the logarithm in ( ^2]) in powers of ^p and s. For ^ = this shift eliminates 
the odd powers of s (of order higher than unity) from the expansion, thereby simplifying 
subsequent calculations. We will truncate the expansion at second order in ip, thus treating 

11 



the electrostatic interactions only at the DH limiting-law level. The expansion in s will be 
retained to fourth order: we assume s(r) to be the basic order-parameter field and, as in the 
usual LGW formulation,^'E^ this is the minimum number of terms required to study simple 
critical points. The grand partition function then reads 






(24) 



*j 



where we have put /i = 1 + 2 (ln2 — /i) ^ . ipj^{0,j). 

In order to perform the integration over the non-critical field ip, we rewrite the partition 
function in reciprocal space, replacing lattice sums by integrals, using 

'^ n k 

where f2 is the first Brillouin zone of the reciprocal lattice.Edfill We further suppose that the 
Fourier transformed short-range potential is of the formE2l 

0s{k)-AJo{l-Rie), (26) 

where Jq is related to the strength of the interaction (divided by ks T) and Rq to the 
interaction range, both of which can be determined for a given model. The values of Jq 
and Rq affect the nonuniversal critical parameters (temperature, chemical potential and 
crossover temperature) but, as we will see below, estimates of the crossover temperature, tx, 
may be obtained by examining various typical values for Rq. As in the discussion following 
(|I^), we have neglected in (^Bj) higher order lattice terms, etc., of the form fc^a^, etc. 

Finally, in terms of the ionicity X, defined in (|T]) it is convenient to introduce the small 
perturbation parameter y and to rescale the ip field in reciprocal space according to 

/ = 4.I.1J, ^W = M. (27) 

[Since we are interested only in the critical region, we may, here, neglect the difference 
between T* and T*, as used in (|l]).] The reduced Hamiltonian, —H/ksT, in (plf) then 
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becomes 

k fci ^2 

~ 12 / / / / ^^^^^ ^^^"^^ '^^^^^ ^^^^^ V^(-^i -h-h- h) - Afisrs{0) (^28) 

fcl fc2 fcs ^=4 

- 2 / ('^sr- + r^^r ^^) \K^)\'^ -'^sr j j j s{ki) s{k2) s{h) s{-ki -k2- /cs), 

k fcl fc2 fcs 

where Afisr = 1 + (ln2 — /i)/2 Jq while 



^-1, rl = ^, «.. = 1 (29) 



are the coefficients of the reduced Hamiltonian of the short-range uncharged fluid. 

Now on integrating over ip in the partition function (0), the coefficients of s" in Ti 
become renormalized because the fields are coupled through the cubic and fifth order terms, 
sip'^ and s^ip"^. This integration over ip as a simple Gaussian field could be carried out 
exactly but it involves a nontrivial diagonalization of 7i which is not required, since the 
resulting Hamiltonian is to be expanded only to fourth order in s, for consistency with the 
previous truncation. (However, the study of BrilliantovEj is instructive in this respect.) In 
fact, the simplest way to calculate the electrostatic corrections to the coefficients of s", is 
via perturbation theory using Feynman diagrams. To this end we write the partial reduced 
Hamiltonian 

,2 „,2 



n^ = -\j{\ + k')\m? + ^-^jjs{h)i^{k,)i;{-k,-k,) 



k fci,fc2 

,2 




T^ 1 1 1 I -^(^i) -^(^2) siks) ^{ki) ip{-ki -k2-k^- k^), 



(30) 



fcl fc2 fca fc4 

where the ffist term is simply Gaussian while we regard the second and third as the "inter- 
action" . 

A standard cumulant expansionE^'Eil yields the n-point correlation functions generated 
by the Gaussian propagator. By carrying the expansion out to fourth order, we obtain a 
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Hainiltonian for a single fluctuating field in the form 

TYo = -- / (ro + Tg A;^) |sp - vq // sss - uq /// ssss - A/iq sq. (31) 

k fcl fc2 fel ^2 fcs 

where the coefficients, renormalized by the electrostatic interactions, are given to leading 
orders by 

ro = rsr- ly^hiy), (32) 

r^ = rl + ly'l3Aiy), (33) 

v, = ±y'h{y)-±y^I,{y), (34) 

Uo = ^[l + ly'l2{y) -|//4(l/)], (35) 

A/io = A/i,,-i|/2j^(^)_ (36) 

If the first Brillouin zone of the reciprocal lattice is approximated by a sphere of radius re /a 
the propagator integrals, In{y), are given by 

TV 



2 7r2 7 (a;2 + 7/2/2)" 





'3,4-2vr2 ] '^''l(x2 + y2/2)3 3 (x^ + |/2/2)4 T ^^^^ 



Clearly, when the ionicity vanishes, y ^ 0, the coefficients r, t, ■ ■ ■ reduce to those of the 
short-range, uncharged fluid. 



Finally, the third order term in ( plj) may be eliminated in the usual way, s =^ s+Vq/A uq, 



yielding a standard short-range LGW Hamiltonian, namely, 

Hs = -- hr + T^ P) \s\^ - u /// ssss-A^iso, (39) 



fci k2 k-i 
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with coefficients 



r = ro-f^, r^ = rl u = u^, A/x = A^o - ^ + ^. (40) 

4 Mo 4mo 8uq 

Accepting the hmitations imposed and the approximations made, we are now in a position 
to study the effects of the electrostatic interactions on the critical behavior of a simple fluid. 
The easiest way to proceed, as indicated above, is to use the Ginzburg criterion to estimate 
the crossover temperature, tx, as a function of the ionicity X oc y"^. This criterion indicates 
the size of the temperature interval around the critical point outside which a mean-field 
description is self-consistent. For temperatures inside that interval the mean-field description 
breaks down and Ising behavior (in this case) should manifest itself. The temperature scale 
so set has become known as the Ginzburg temperature, tc, and it is reasonable to take it as 
an estimate of the crossover temperature, tx ,t3T^ although the presence of significant higher- 
order terms in (|28| ) could undermine the practical validity of the criterion.li^cfl In the next 
section we compare the Ginzburg temperatures of systems with finite ionicity (X oc y^ > 0) 
and without electrostatic interactions {y = 0). 

IV. APPLICATION OF THE GINZBURG CRITERION 

The condition for a self-consistent mean-field description is that fluctuations of the field 
(calculated using the mean field) are much smaller than the typical or average mean field. 
For the s^ Hamiltonian (^) studied below Tc, this condition may be expressed conveniently 
all 

3^2 u ^ ,,^, 

Naturally, in deriving this, we have set AyU = 0, so that the critical point is approached 
along an axis of asymptotic symmetry. It is then appropriate to note that the short-range 
coefficient, r^^, can normally be written asE^I 

-'-cfl 
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where TJ^q is the mean-field critical temperature of the uncharged system. For a system 
with electrostatic interactions the mean-field critical temperature, Tcfi{y), increases by an 
amount which follows from (^2|) and (|40| ) as 



tM^If^-l = ly*IM + ^J^. (43) 



'c,0 

We may use this result to rewrite the coefficient r in (^) as 

T - T,,o(y) 



T'o 



l + to{y)]t, (44) 



where t = {T — Tcfi)/Tcfi is the (mean-field) reduced temperature of the charged system. 
The Ginzburg criterion then becomes 

,_, 18 u^iy) 

"^^ ^ 7r2 [l + to(y)]r6(|/) ^ ^ 

which defines the reduced Ginzburg temperature, to- At this point the value of to depends 
on the ionicity X = y^/47r through T{y), u{y) and to{y), as given in (P5D, (P5|), ( |iD| ) and (^3|), 
respectively, and on the range, Rq through Tgr- The dependence of to on Jq has dropped 
out since, using (P^D, one finds r^^^ = [1 + toiv)] Rq- 

For the short-range uncharged system the Ginzburg temperature reduces to 

If, as a convenient and not unreasonable value, one adopts Ro/a ^ 0.4828, one obtains 
tci^) = 1- For this particular short-range model, we have calculated the parameters of the 
reduced Hamiltionian (|39D of the system with electrostatic interactions, as a function of the 
ionicity, I. The results for r, u and the shift in the mean field critical temperature of the 
system are plotted in Fig. 1, as functions of I. The corresponding Ginzburg temperature, 
tci^), is displayed in Fig. 2. The effective interaction range, r increases monotonically with 
ionicity while the coefficient, u, of the fourth order term exhibits a maximum at X ~ 3 
and decreases slowly at higher ionicities. However, within our Debye-Hiickel limiting-law 
approximation, u remains positive for all ionicities with, in fact, u ^ n(l + n^) — 0.105 as 
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y -^ oo. As a result the Ginzburg temperature, tc-, varies nonmonotonically with ionicity, 
1: see Fig. 2. A maximum, to — 1.03, occurs at X ~ 0.7 indicating that, in general, the 
electrostatic interactions reduce (by up to ':^ 30%) the interval around Tc in which Ising 
critical exponents may be observed. Closer inspection of the behavior of tc with I reveals 
an initial drop off of tc-, as shown in the inset. 

In order to check the robustness of these results, we have calculated the ratio of reduced 
Ginzburg temperatures, tG{y)/tG{0), for other values of Ro/a in the range 0.1 to 10 (cor- 
responding to uncharged Ginzburg temperatures, ^^(O), spanning 12 orders of magnitude). 
The results are plotted in Fig. 3. The behavior of the ratio of Ginzburg temperatures, 
tG'(y)/i^G(0), varies with the particular value of Rq but at ionicities X > 2, a reduction of to 
is found for all the systems. At low ionicities, however, the behavior of to depends on the 
range of the nonionic interactions. The initial drop off is more rapid for low values of Rq/q. 
For the system with R^/a = 0.1 a minimum value of to is clearly visible. 

In assessing the results shown in Figs. 1-3, the perturbative nature of our calculations 
must be borne in mind. The scale of probable validity in X may be guaged by noting that 
Monte Carlo results for the RPM,E3E3iHZl the extreme limit of "complete ionicity," indicate 
X ~ 20.4. Clearly our results are, at best, of qualitative validity when X exceeds, say, 5. 

Finally, some contact with experimental results may be made by estimating the values 
of X using the observed reduced critical temperatures and other data: see Table || which 
presents the Ginzburg temperatures calculated for various experimental systems (assuming, 
for simplicity, R^/a = 0.4828) and the corresponding, experimentally assessecS crossover 
temperatures, t^- Realistic values of X lead to the prediction of reductions of the reduced 
Ginzburg temperature that correlate with the experimental data. However, our theoretical 
results for to grossly underestimate the reductions in tx observed in the real ionic systems 
listed. 
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V. DISCUSSION AND CONCLUSIONS 

Starting with a Hamiltonian with attractive short-range as well as long-range Coulomb 
interactions, we have derived a field-theoretic description for a symmetric model with two 
coupled fiuctuating fields, ip (conjugate to the charge density) and s (conjugate to the 
number density). This lays a foundation on which to study explicitly the infiuence on 
critical properties of the ionicity, that is the strength of the Coulomb coupling relative to 
the short-range interactions. In an exploratory approach we expanded the weight function 
to second order in ip, corresponding to a Debye-Hiickel limiting-law level, and to fourth order 
in s, the usual LGW level. Finally, we integrated out the ip field so generating electrostatic 
corrections to the coefficients of s" in the final effective Hamiltonian, Tig — see (^) — in a 
systematic way. 

We examined the perturbative effects of the electrostatic interactions on the near-critical 
behavior of a simple fluid by using the Ginzburg criterion to estimate the crossover temper- 
ature, tx, as a function of the ionicity, deflned more explicitly by X = e^/D a ks T: see (p. 
The calculations indicate that as X increases, the effective range of interaction, r, increases 
slowly while the coupling of the fourth order term in Tig-, namely m, decreases for realistic val- 
ues of X and Rq after an initial, relatively rapid transient. As a result, the relative Ginzburg 
scale, tG(2^)/^G(0), [see (^)] decreases with increasing X, by up to 20 — 30% (within the 
approximations): see Fig. 3. This suggests that the expected domain of Ising-type behavior 
around Tc should, likewise, be reduced by the presence of Coulombic forces. 

Although we believe these results are interesting, one must certainly be cautious in 
applying them to real ionic systems. First, it is not certain that the Ginzburg criterion 
provides a sufficiently reliable measure of tx-l^rO£ltJ" Indeed, despite the correct trend of to 
with the ionicity X seen in Table II, the predicted magnitudes of the reductions of tx are 
vastly smaller than observed in the real systems studied. 

Second, in our perturbative treatment of the model the critical behavior is always driven 
by the short-range interactions, and only nonuniversal critical parameters are affected by 



the electrostatic interactions. It is unlikely that such a mechanism describes satisfactorily 
the criticality of all real ionic systems. In fact, as argued by Pitzerjl^ there are chemical 
systems in which the ionic interactions are almost surely the principal determinants of the 
critical behavior. And, of course, in the restricted primitive model (RPM) only Coulombic 
forces act but a critical point is surely realized.Q'tHcIl Such systems clearly lie outside the 
range of validity of our truncated, perturbative treatment. Furthermore, as suggested both 
theoreticallyEfE^ and by the experiments of Pitzer and Narayanan,ll3 Coulombic electrolytes 
of high ionicity may lie on a different global thermodynamic locus than do solvophobic 
systems such as we have sought to describe in our present approach. 

Within the approximation presented, u never vanishes; however, it is possible that carry- 
ing the calculations to higher ordercJ would lead to negative u for sufficiently large ionicity. 
At a mean-field level in a simple lattice system that would suggest tricriticalityu and the 
presence of an "antiferromagnetic" or crystal phase. But, by any measure, one would then 
need to retain sixth and, possibly, higher order terms in the effective Hamiltonian;E^ and, in 
that case, a renormalization group treatment would seem necessary to elucidate the nature 
of the resulting transition(s).c3'Ea It should be forcibly stressed, however ,0 that a tricritical 
point per se cannot arise in a normal one- or two-component fluid system except by some 
special, nongeneric accident! (In a three-component system a tricritical point can occur but 
is characterized by a particular temperature, Tj, and three flxed densities, say, p^, p^ and 

The present approach has been extendedE^I to asymmetric models with distinct short- 
range interactions, <f''^^_^{r), (f'^^_{r) and (f'L_{r). The calculations yield an effective Hamil- 
tonian very similar to that obtained by Nelson and FisherEa in their study of metamagnets 
in the presence of an external magnetic fleld (which, when large enough, induces tricritical- 
ity). A renormalization group calculation employing the e-expansion has been performed:EHl 
it yields an unstable Gaussian flxed point, a symmetric Ising-type (Fisher- Wilson ?T,=l)EyJ^ 
flxed point, and two further asymmetric Ising flxed points. In addition, a separatrix appears, 
at the s^ level of approximation, that delimits a region of "runaway" flows: such behavior 
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requires proper study of the Hamiltonian with higher order terms retained.E3 However, as 
hinted above, the significance of such flows within the context of fluid systems is unclear at 
this stage. 
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TABLES 
TABLE L Experimental melting and critical temperatures, T^ and Tc, in various solvents for 
three organic salts: (a) tetra-n-pentylammonium bromide; (b) triethyl n-hexyl ammonium triethyl 
n-hexyl boride; (c) tetra-n-butylammonium picrate. 

Salt Tm (K) Solvent T^ (K) 

(a) ~ 374.05 waterS 404.90 

(b) ~ 251 diphenyl etheiili 309-312 

(c) 364-365 1-choloroheptanei 404.4 

1-tridecanolil 342 

l-dodecanolO ~ 332 

1,4-butanedioiii ~ 333 



TABLE IL Experimentally assessed crossover temperatures, tx, and reduced Ginzburg tem- 
peratures, tc, calculated assuming an uncharged model with Ro/a = 0.4828 [or tciO) = 1] foic2l (a) 
tetra-n-butylammonium picrate (TBAP) in l-tridecanol;Ej (b) TBAP in l-dodecanol;Ej (c) TBAP 
in 75% 1-dodecanol plus 25% l,4-butanediol;Ej (d) Na in NH3;ll3a (e) tetra-n-pentylammonium 
bromide in water J13 

System ionicity, I tx to 

uncharged fluid 0(1) 1 

(a) 17.9 ~ IQ-^ ~ 0.712 

(b) 16.8 ~ 0.9 X 10-2 ~ 0.717 

(c) 8.9 ~ 3 X 10-2 r^ 0.777 

(d) 6.97 0.6 X 10-2 ^ 0.807 

(e) ~ 1.4 0(1) ~ 1 
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FIGURES 
FIG. 1. (a) The effective interaction range, r, in the Hamiltonian ( |39| ) as a function of the 

ionicity, X, for a system with R^/a = 0.4828; (b) the fourth order couphng constant, u, in (^) as 

a function of the ionicity; (c) the reduced shift of the mean-field critical temperature, to- Notice 

that both u and to are independent of the value of Ro/a. 

FIG. 2. The reduced Ginzburg temperature to as a function of ionicity, I, for the short-range 
model of Fig. 1, which is described by tciO) = 1. The inset shows the behavior of to close to the 
origin. 

FIG. 3. Ratio of the Ginzburg temperatures calculated for a system of given ionicity to those 
for the corresponding short-range models tG{'I)/tG{0) for the cases Ro/a = 0.1 (dashed line), = 0.2 
(dashed-dotted hue), = 0.4828 (full line), and = 10 (dotted line). 
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